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We show that the Stephenson-Kihuister-Yang (SKY) equation combined with Camenzind's matter 
current term naturally provides the cosmological constant and dark radiation as integration con- 
stants of the SKYC field equation. To characterize the property of the dark radiation, we develop 
a method to separate it from the ordinary radiation. We found a special property of Camenzind's 
matter current, namely that the solution space for radiation in fact belongs to that of the vacuum 
solution of SKY equation. We also found that his matter current does not obey the conservation 
condition suggested by Kilmister. Finally, we discuss the possible role of dark radiation emergent 
from the SKYC theory in recent cosmic-microwave-background observations and its implications to 
the inflation scenario. 
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I. INTRODUCTION 

Various attempts have been made throughout the last 
century to unify all fundamental interactions. One pio- 
neering effort was made by Weyl in 1918 [ll, [3> where he 
assumed that all physical laws should be invariant under 
conformal transformation. This seminal viewpoint intro- 
duced by Weyl is now known as the gauge invariance. 
Based on the principle of gauge invariance, Weyl refor- 
mulated the connection in Riemannian geometry to unify 
gravity and electromagnetism. 

The gauge theory of Yang and Mills has inspired a new 
type of relationship between geometry and physics [3|. 
Such a theory can be considered as a vector bundle on a 
Riemannian manifold [^-l6|, where a section of a vector 
bundle corresponds to a matter field and a connection on 
the Riemannian manifold corresponds to a gauge field. A 
gauge theory can therefore be recognized as a functional 
(called Yang-Mills functional) that acts on a metric con- 
nection on a vector bundle. Yang-Mills functional is also 
invariant under gauge transformation. Thus, Yang-Mills 
theory exhibits a close correspondence with the vector 
bundle theory in differential geometry. 

We now know that all fundamental interactions: the 
electromagnetic, the weak, and the strong, can be de- 
scribed in the language of gauge theory except general 
relativity (GR). The pioneer works to formulate Ein- 
stein's GR into a gauge theoretical framework started 



* |pisinchen@phys.ntu.ed u.tw] 
t izumi@phys.ntu.edu.tw 
^ ■d96222003@ntu.edu.tw. 



with Utiyama, who suggested that GR can be written in 
the language of the gauge theory if the symmetry group 
is chosen either as the Poincare group or the translational 
gauge group [THlOj . 

Yang proposed his gauge theory of gravity based on 
the GL{n) group and identified the affine connection as 
the gauge potential ^llj], with the action for the pure 
space proportional to R'^'^'^^R^^aX, where R^^aX is the 
Riemann tensor. Prior to Yang's proposal, Stephenson, 
Kilmister and Newman had derived a similar equation 
but with additional constraints from a second-ordered 
curvature Lagrange density [IJ, [lj| . We shall call their 
field equation the SKY equation. The SKY equation 
reproduces all solutions of the vacuum Einstein equa- 
tion |12h17I|. Later Camenzind proposed a matter cur- 
rent term for the SKY field equation [1^, although it 
was not deduced from an action. This task was fulfilled 
by Cook in 2009 [19]. We shall call the complete theory 
that includes both pure space and matter contributions 
the SKYC gravity. A renormalizable quantum theory of 
gravity based on the Einstein-Hilbert action is known to 
be difficult to attain |20l - t22l |. Being a gauge theory, a 
quantum gravity based on the SKYC theory mig ht stand 
a better chance to be renormalizable [111, [2J, [2J| . On the 
other hand, the high-order, quadratic contraction of Rie- 
mann tensor without the protection of the Gauss-Bonnet 
condition will necessarily render the ghost problem. One 
may, however, work around the ghost problem by ex- 
pelling the negative pole to the UV limit, say to the 
Planck scale, and regarding the SKYC gravity as an ef- 
fective theory |25| . 

There has been the longstanding cosmological constant 
problem [2a|) The nature of the cosmological constant 
(CC) term introduced by Einstein to his field equation 



is a priori undefined. On the other hand the quantum 
vacuum energy satisfies the properties of the CC, yet its 
value is about 124 orders of magnitude larger than the 
critical density of the universe, which is comparable to 
what is required for CC to explain the observed acceler- 
ating expansion of the universe |27l . |28| . From the math- 
ematical consideration, the CC term cannot be removed 
from the Einstein-Hilbert action. The SKYC equation is 
second order in terms of the afEne connection. While the 
metric is a priori not a dynamical variable, in order to 
reduce the SKYC equation to the Einstein equation one 
define the relation between the connection and the metric 
in the usual way. As a result the SKYC equation can be 
redressed as a 3rd order differential equation of the metric 
where the CC term is necessarily absent. The CC term, 
however, is recovered as an integration constant in re- 
ducing the SKYC equation to Einstein equation. In this 
approach CC is no longer arbitrary but determined by 
the boundary condition of the universe, which is geomet- 
rical in nature and has nothing to do with the quantum 
vacuum energy. For example it has been proposed that 
the underlying geometry is the de Sitter space and this 
integration constant has been associated with the radius 
of curvature of the de Sitter space [13, ^M, 13 ■ Therefore 
the SKYC formulation of gravity may provide a solution 
to the CC problem. It goes without saying, of course, 
that a big part of the problem still remains: why does 
the quantum vacuum energy not contribute to gravita- 
tion in the first place? This issue is beyond the scope of 
this paper and we will not dwell on it further. 

Aside from the CC problem, recent observation data 
indicate that there is an apparent increase of the effect 
number of neutrino flavors, N^s, between the big bang 
nucleothynthesis (BBN)and the cosmic microwave back- 
ground (CMB) epochs [30|-[3J|. The standard model pre- 
dicts Ncs — 3.04 at the epoch where the universe was 
dominated by photons and neutrinos after the electron- 
positron annihilation. However, recent experimental 
data suggest a larger A^eff |30l - l34l | . It has been suggested 
that the difference, AA'cg = A'eff — 3 ~ 1 at nearly 2a 
level, may indicate the existence of an additional, hereto- 
fore unobserved "dark radiation (DR)" density. On theo- 
retical side, this additional radiation was also called upon 
by the brane- world inspired and other cosmological mod- 
els |35l - l4l| . It happens that another integration constant 
in the SKYC theory has exactly the required character of 
dark radiation. We note, however, as has been pointed 
out recently by Birrell et al. [42|, that the fact that 
neutrinos have rest mass and that their distributions are 
non-thermal under free-steaming, can well explain such 
an increase of N^s without the need to invoke dark radi- 
ation. With that in mind, the dark radiation arisen from 
SKYC gravity is nonetheless a free bonus at our disposal 
subject to observational constraints. 

The purpose of this paper is to examine explicitly the 
salient features, in particular the CC and dark radiation, 
of the SKYC gravity mentioned above. To do so, we de- 
rive the Friedmann-Lemaitre- Robertson- Walker (FLRW) 



equation associated with the SKYC gravity with integra- 
tion constants. We identify one integration constant as 
CC and the other as dark radiation. Subsequently, we 
develop a method to investigate the property of dark ra- 
diation and demonstrate that the density of dark radia- 
tion is a constant on the constant-time hypersurface in 
every metric. 

This paper is arranged as the following. In Sec. |IT1 we 
review the SKYC gravity. In the next section, the SKYC 
field equation on the FLRW metric is derived and inves- 
tigated. We discuss the property of nullity in Camen- 
zind's matter field in the perfect fluid model. In Sec. IIV[ 
we characterize the properties of the dark radiation. In 
Sec. |Vl we summarize and discuss our findings. In the 
appendix, we discuss the problem with Cook's theory, 
which intends to give rise to the SKYC equation from 
the action level. 



II. SKYC EQUATION 

In a non-abelian gauge theory, the field strength is de- 
fined as 



F„, 



-[^M.-D.], 



(1) 



where I?^ = d^—igA^ and A^ is the gauge potential. The 
field strength F^^, can therefore be expressed in terms of 
the gauge potential Ay as 



df,Ay - d„A^, ~ ig[A^,,Ay]. 



(2) 



When connecting Riemmanian geometry with the non- 
abehan gauge theory, the Christoffel symbol, T^^, and 
the Riemann tensor, R^va\, play the role of the gauge 
potential ^^ and the field strength F^,y , respectively. Mo- 
tivated by this gauge connection, Stephenson, Kilmister 
and Newman investigated an alternative theory of grav- 
ity and obtained the following equation of motion. 



v^Rr 



fj.al3 



0, 



(3) 



which is a hig her-order derivative with respect to the 
metric [l^, [l^] . Invoking the second Bianchi identity, it 
can be readily verified that Eq.® is equivalent to 



V^i?a^ — V^i?Q^ — 0. 



(4) 



In fact, Eq.([3l) is exactly Yang's gravitational field equa- 
tion for pure space under the GL{n) gauge symmetry 

%■ 

To complete the gravitational field equation, Kilmister 
introduced the matter current / such that |13l] 



^aRpfj. — ^fsRafi — lafin- 



(5) 



He found that / has 20 independent components under 
the constraint of the conservation law: 



y^ja/3M = 0. 



(6) 



However, he did not provide the exphcit form for the 
current density I . 

Later, Camenzind wrote down a Yang-Mills field equa- 
tion for 5*0(3, 1) with a current J different from Kilmis- 
ter's /: 



where the current density J"*^'' has the form: 



(7) 



(8) 



V" (T^'^ - l/2g''^T 



Here T^" is the energy-momentum tensor [l8|. Never- 
theless, an action associated with this matter current re- 
mained lacking. Aspired by the analogy with Maxwell's 
theory. Cook later proposed an action term for the mat- 
ter current, F^oJ"'', from which the SKYC field equa- 
tion can be derived [I^ • It can be verified that solutions 
of Eq. © cover the entire solution-space of the Einstein 
equation with source |43|. 



III. SKYC EQUATION IN FLRW METRIC 

In this section we will use the FLRW metric to de- 
rive a modified Friedmann equation from the SKYC field 
equation, Eq.©, and discuss the nullity of Camezind's 
current density in the radiation-dominated case. 



A. EOM in Homogeneous Universe 

We consider a homogeneous and isotropic universe and 
use the FLRW metric to study Eq.©. The metric is: 



where no summation over the repeated index "z" occurs. 
Other components vanish due to the symmetry of the 
metric. The expression for T^i,, Ea. (fTO|l . leads to 

a^'a + aha — 2d^ — 2ka = SttG ( tto (p — p) + a{p + p) ) , 



which can also be expressed as 



(14) 



H + mH -2k^ = ^ttG C^{p - p) + H{p + p) ) . (15) 



Eg. ([15)) is the modified Friedmann equation in SKYC 
theory. By integrating Eq. p^ with Eq. p^ . we can ar- 
rive at 



H + 2H^ + ^ = 



SttG /I 



gP-P 



2Gi, 



(16) 



where Ci is an integration constant. By performing an- 
other integration to Ea. (fT6| . one obtains 



H' = 



SttG 



k ^ G2 

— + Gi + — . 



(17) 



This is the Friedmann equation for the SKYC theory. 

Eq. pT)) can be interpreted as follows. In a FLRW uni- 
verse, Gi plays the role of the CC term. As will be shown 
in the following section, the constant Gi plays the same 
role as the CC even in an inhomogeneous universe; this 
fact has been shown by Cook [19]. (Our Gi is identical 
to the trace part of A^i^ in his article.) The term involv- 
ing G2, on the other hand, scales the same as radiation. 
Since it does not couple with any matter except grav- 
ity, it can be identified as the dark radiation. We note 
that although the G2 term is not induced by the matter 
current, it does affect the evolution of the universe. 



ds^ = -dr +a^{t) 



dr^ 



1 — fcr^ 



r'^dO?- 



(9) 



where a{t) is the scale factor. From the symmetry prop- 
erty of the FLRW metric, the energy-momentum tensor 
T^j/ can be written as 



Tp,„ = {p + p)u^u^ +pgf^u, 



(10) 



where u^ = {d/dt)^. In Eq.©, we assume T^j, to satisfy 
the conservation law: 



V^Tf"" = 0, 



(11) 



which in turn leads to the usual equation 

p^-3H{p + p), (12) 

where H = a/ a. 

The only nontrivial components of Eq.© are: 



^tj^^ioi — ^TrJoii, 



(13) 



B. Nullity in Current Density Tensor 

In this subsection we show that Camenzind's current 
density will be null if the equation of state for radiation 
is invoked. On the other hand, such equation of state 
belongs to the vacuum equation in the SKYC gravity. 
We study this special property here. 

The equation of state for radiation is p = (l/3)p. 
Based on this. Eg. ([12)) becomes 



Prtt 



E. 



where E is & constant. Thus the dynamical equation for 
a is 



SttGE 
3a2 



-k + Cxc? 



9i 

a2' 



(18) 



where the first term on the RHS represents the contribu- 
tion from the ordinary radiation. As we have commented 
in the previous subsection, the Ci term also behaves like 
radiation. Hence for a perfectly homogeneous universe 



one cannot distinguish the ordinary radiation from the 
dark radiation. 

There is a property about the nullity of the current 
density tensor to be mentioned. If one inserts Eq. p^ 
into the gravitational equation Eq.©, then the RHS of 
Eq.© becomes null, that is, 



^uR naB — 0, 



(19) 



yet the current density is a priori not necessarily zero. 
To avoid the contradiction, one must impose the ansatz 
that Japii equals to zero in the radiation case. Under the 
assumption of p = wp, we find that both w = —1 and 1/3 
will make JaPfi, vanish. In this situation, Eq. (1181) now be- 
longs to the set of the vacuum solutions of Eq.©. Such a 
higher-order gravitational equation brings the radiation 
terms into integration constants of the vacuum equation 
of Eq.® from the current density Jap^- This is an in- 
teresting property in Camenzind's current density of the 
matter. 

Although the genesis of the dark components from 
Eq.© is interesting, it might be a signal of an incon- 
sistency. This property means that Camenzind's current 
density for radiation and dark energy cannot be a source 
of gravity in the exact FLRW universe. This result stems 
from the assumption in Eq.®, where the energy density 
T^i, is introduced, that the energy density T^i, gives rise 
to the current and that this current satisfies the conser- 
vation law. A deeper inspection shows, however, that 
while the Camenzind's current density J%^ is the cur- 
rent for gravity, yet there is no theoretical origin for the 
energy conservation law. Furthermore, the non-trivial as- 
sumption of Eq.® might result in the emergence of the 
dark components and the fictitious property about nul- 
lity because Camenzind's current density is in the form 
of higher order derivatives, although this assumption is 
needed to recover GR. 



IV. CHARACTERISTIC OF DARK RADIATION 

In the previous section, we have seen the fact that one 
can not separate the effect of the dark radiation from that 
of the ordinary radiation in the homogeneous universe. 
In order to distinguish them, we consider the case of an 
inhomogeneous universe. 

Firstly, we solve the EOM in the inhomogeneous uni- 
verse. We define the effective energy-momentum tensor 
as 



Because of Eq. ([TT|) , the divergence of this gives 



V^T^, = 0. 



(20) 



(21) 



Substituting ([?T|) into the gravitational EOM, we have 
VaT^^ - V^Tq.^ - - [gp^Vaf - ga^V pfj = 0,(22) 



where T is the trace of T^^ . Multiplying g'^'^ by Eq. 
we can obtain 



^pT = 0, (23) 

where we use Ea. (l?I|) . Integrating this, we have 

f = 4Ae//, (24) 

where Ae// is a constant. As we have commented at the 
end of Sec. IIII Al it plays exactly the role of CC. 

We now separate the contributions other than the ef- 
fective CC from T^^: 



^ iiu — ^ ^u ^^cffg^u- 



(25) 



We know from Eqs. (I?T|) and (|^^ that the following equa- 
tions must be satisfied: 



5 = 0, 


(26) 


v^s^, = 0, 


(27) 


Va5;37 - V/J^o^ = 0. 


(28) 



If the form of S^^ is the same as that of radiation fluid on 
FLRW metric, then it can be shown that these equations 
are satisfied. Therefore, 3^^,^ is related to the C2 term. 

In order to see the difference of the effective energy- 
momentum tensor S^u from a real radiation fluid, we 
analyze its property on general metric in the form of a 
fluid: 



Sfiu = (p + p)Uf^Uy + pgf^u, 



with 



1. 



(29) 



(30) 



The traceless condition (l^5|) fixes the relation between p 
and p as 



p-\p- 



Then S^n becomes 



where 



(31) 

(32) 
(33) 

(34) 

is the induced metric on the hypersurface which is or- 
thogonal to u^. 

Eq. ([77]) can be written as 

4 14 

-upuf'd^p + -dfip + -p {Uf.Wuf} + upWuf,) = 0. (35) 

Multiplying u^ and hi^ by the above expression, respec- 
tively, we can obtain 

4 



Q 4 1 

= pUf^u^ + -phfj^^, 



'^fiL^ 9fiL^ ~r ^/j.^t/ 



u^'d^,p = --pV^u^, 
h/dfsp = -ApUf^Vu^, 



(36) 

(37) 



where we use 



/V^u^ = iv (u/^up) = 0, 



(38) 



hJVufj ^ (Sj^ + u^u^] Vufj = Vu^. (39) 

Multiplying u^ by Eq. ([5S)) and subtracting it from 
Ea.([571). we have 






(40) 



On the other hand, Eg. ([28]) can be written as 

UpU^daP - UaU^dpp + -hp^daP - -ha^dpp 

4 

— M-yVcU^ + U-yV^Uct) — 0. (41) 

Multiplying u'u^ by Ea. (PT|) . we have 



hjdfip - -pu^XI pUa = 0. 



(42) 

Combining it with Ea. (|37l) . we have 

u^V^u, = 0, (43) 

hjdup = 0. (44) 

Multiplying hy^u^ by Eq. (|lT|) . we have 

- -^XaU^dfiP - X/OVaWA = 0, (45) 

where we have used Eq. (j^^]) . Combining it with Eq. (IMl) , 
we find 



VaWA = -ha\V^Uf_,. 



(46) 



Multiplying u^h\^ and hf^h\^ by Eq. (PT|) . respectively, 
gives 



(47) 



Uahx'^di3P+ -pVaUx 

--p'^XUa - -Uxul^V/SUa = 0, 
-h^xdaP - -heahx^dpp - -pUaVxUe = 0. (48) 

Substituing Eqs.gO]), (031), gH) and (ge]) into Eas. (l47|) 
and (H5)) . we see that they are automatically satisfied 
and no additional condition is obtained. 

In summary, we have transcribed the original equations 
for Sf,^ into Eqs.(l36l), (p)) . (p ]) . (|46l). and 



9wP = ■rpWi^V'u^i 



(49) 



that govern the characters of p and w^. In turn, these 
conditions constrain the form of the effective energy- 
momentum tensor. Equation (j37p means that energy 
density must be constant on the hypersurface which is 
orthogonal to u^. Therefore, the dark radiation from 
the effective energy-momentum tensor can affect only the 
background dynamics and we confirm the existence of the 
real radiation fluid from its perturbation. 



V. CONCLUSION 

We investigated the SKYC theory of gravity by way 
of solving its field equation, Eq. d?]) , in a FLRW universe 
and arrived at a modified SKYC Friedmann equation, 
Ea. (|15p . Being a higher order derivative equation than 
that for GR, the SKYC Firedmann equation gives rise to 
two integration constants when it is reduced to lower or- 
der. One of the two is clearly related to the cosmological 
constant while the other is related to the dark radiation. 
We have also demonstrated, in a homogeneous universe, 
that this dark radiation is indistinguishable from the or- 
dinary radiation. In addition, we pointed out the nullity 
of the current density J in the radiation case under the 
FLRW metric. 

In order to further pin down the nature of our dark ra- 
diation, we turned to a general, inhomogeneous universe 
and introduced a methodology to look for its possible dif- 
ference from the ordinary radiation. We solved the EOM 
in the inhomogeneous universe and constrained the form 
of the tensor S'pi, in Sec. IIVI We found that if p is fixed at 
one point, then it will be a constant on the hypersurface 
orthogonal to u'^. That means that there does not exist 
any degree of freedom for the perturbed dark radiation. 
That is, this SKYC dark radiation only has the zero mode 
term but has no perturbed term. In contrast, the ordi- 
nary radiation can be perturbed and can therefore prop- 
agate in all spacetime. We conclude that the SKYC dark 
radiation is indeed different from the ordinary one. We 
should like to comment, however, that in our derivation 
in Section HVl the expression 5*^^ = [p + p)u^Ui, +pg^i, 
is not the most general form. We will pursue a more 
general expression for it in our future work. 

Some comments are in order with regard to the rela- 
tionship between the SKYC dark radiation and the infla- 
tion. The dark radiation C2 term is an integration con- 
stant in this theory. That is, it is determined by the ini- 
tial or boundary condition of the universe. If the SKYC 
theory is incorporated with inflation, then the density 
of the dark radiation pdr must start from a tiny value. 
Otherwise the inflation can not be trigged because pdr 
scales as a~^ and dominates at early times. At the end 
of inflation and after ~ 60 e- foldings, the scale factor has 
grown by ~ 10^'^ times. This means that pdr must be 
smaller by 80 orders of magnitude than poR at late times. 
Being so tiny, we may as well set C2 to zero. On the other 
hand, if SKYC theory does not include the inflation sce- 
nario, then the C2 term can in principle be identified with 
the dark radiation and be fixed by the observation data. 
The recent Planck data gives TVoff = 3.36+j]j|(95%) 
based on the combination of WMAP -|- highL data [3J] . 
However, this fit produces a 2.5 s.d. tension with di- 
rect astrophysical measurements of the Hubble constant. 
Including priors from SN surveys removes this tension 
and results in iVeff = 3.62l°;^g(95%). The larger iVeff 
suggests a need for the dark radiation. Without fusing 
SKYC gravity with infiation, we fix C2 with the Planck 
data and find C2 ^ 5.59 x 10"'^^ kg/m'^, while including 



SN priors Planck data gives C2 ^ 1-02 x 10 ^^ kg/m^ 



VI. ACKNOWLEDGEMENT 

It is a pleasure to thank R. J. Adler, Je-An Gu, De- 
baprasad Maity, Yen Chin Ong, Shu-Heng Shao and Yen- 
Wei Liu for helpful discussions. Pisin Chen is supported 
by Taiwan National Science Council under Project No. 
NSC 97-2112-M-002-026-MY3 and by US Department 
of Energy under Contract No. DE-AC03-76SF00515. 
Keisuke Izumi is supported by Taiwan National Science 
Council (TNSC) under Project No. NSC101-2811-M-002- 
103. 



Appendix A: Difference between Stephenson's and 
Yang's Approaches 

Stephenson's gravitational equation is similar to 
Yang's, but there are some important differences between 
the two theories. In brief, Stephenson's theory covers less 
solutions than Yang's even if the connection in Stephen- 
son's theory is identified as the Christoffel symbol. We 
will discuss this point in more details in this appendix. 

In the Palatini formalism, the affine connection and the 
metric are treated as independent variables and the de- 
rived EOMs can determine the relation between the two. 
For example, the relation between the connection and the 
metric can be identified as the condition for the metric 
compatibility when applying the Palatini formalism to 
derive EOMs from the Einstein-Hilbert action. Stephen- 
son and Cook obtained EOMs by applying the Palatini 
formalism to the quadratic curvature Lagrange density 
for pure gravity without matter. About the matter field 
in the SKYC theory. Cook introduced his current density 
tensor at the action level. However, there are problems 
with his approach, about which we will comment in Ap- 
pendix B. 

In Yang's theory, the gravitational force is described by 
the GL(4) gauge field b't. Here the Latin letters stands 
for the indices of the GL{4) gauge group. The action is 



Sb 



dx^V~9 {g^'-g-^ctcu^jip) 



where /° is the field strength: 



fa _ in la _ f~<a ih ic 



(Al) 



(A2) 



and C^^ is the structure constant. 



In order to connect his gauge field b^ to the metric, 
Yang introduced a higher order curvature term in the 
action. The final form of his gravitational action is 



SM'9 



liv\ 



(A3) 



Performing variations of the action with respect to the 
gauge field 6J^ and the metric, one arrives at two types of 
equation of motion. The index a in ¥t has 4x4 values and 
it can be redefined as a = {kl} where k and I run from 

to 3. Yang used an ansatz ha Sk^iSii, = {(^^} that 
satisfies the equation derived from the variation of the 
action with respect to the metric. It gives the relation; 



fiafi) __ pa 

and the other equations become 



(A4) 



(A5) 



This is Yang's gravitational equation based on the gauge 
theory and it covers the solutions of Einstein equation 
for pure space, that is, without matter |ll| . 

In Stephenson's theory, there are two EOMs which 
stem from the variations of the action with respect to 
the metric and the connection, respectively [l2|: 



^Ma/37p 



vafS-y 






M 



2R 



aPui 1 



RaPi^f ~ T.g^R"' RaPfS — 0, (A6) 



Va(i?r% 



0. 



(A7) 



The relation between the metric and the connection, how- 
ever, is different from that in GR. In particular, the affine 
connection in his theory can in principle be different from 
the Levi-Civita connection. If one identifies the affine 
connection as the Levi-Civita connection in Stephenson's 
theory, the resulting EOMs are not equivalent to Yang's. 
Then, Stephenson's equations become 



^ ^uaf3^ 

V,i?7" = 0. 



1 



g^R"^'<'Rap-,s = 0, 



(A8) 
(A9) 



These equations can be recognized as Yang's gravita- 
tional equation fEa. ljASp ) under the additional constraint 
of Eq. (|A9p . In this sense, the solution space of the resul- 
tant Stephenson's equations must be a subset of Yang's 
gravitational equation. Without specifying the connec- 
tion as Levi-Civita connection, Stephenson's equations 
should in principle have different solutions from Yang's. 



Appendix B: Problems with Cook's Theory 

Cook introduced the matter action to Stephenson's 
theory. His recipe, however, not only retains the orig- 
inal problem of Stephenson's theory but also introduces 
another one. The action he proposed is 



Sn^ 



-1 
1677 



(i?" ^"RaPfii^ 



+ 167r J, 



a/3pM 



-g d*X, 



(Bl) 



where the tensor J "^ is Cook's current density, which 
has the same form as Camenzind's. This matter action 
term, however, is not general covariant because the con- 
nection is not a covariant tensor j44| . 

The lack of general covariance must be closely related 
to the non-conservation of the current density J . Ac- 
cording to Noether's theorem, symmetry property of the 



action always goes hand-in-hand with the conservation 
of the current. Therefore, it seems impossible to con- 
struct a general covariant action based on Camenzind's 
current density. This suggests that one should search for 
a different form of the current density other than that 
of Camenzind's, under the constraint that GR must be 
recovered. We will investigate this further. 
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